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Abstract. Characteristic scale lengths of nonthermal X-rays from the SN 1006 NE rim, which are observed by 
Chandra, axe interpreted in the context of the diffusive shock acceleration on the assumption that the observed 
spatial profile of nonthermal X-rays corresponds to that of accelerated electrons with energies of a few tens of 
TeV. To explain the observed scale lengths, we construct two simple models with a test particle approximation, 
where the maximum energy of accelerated electrons is determined by the age of SN 1006 (age-limited model) 
or the energy loss (energy loss-limited model), and constrain the magnetic field configuration and the diffusion 
coefficients of accelerated electrons. When the magnetic field is nearly parallel to the shock normal, the magnetic 
field should be in the range of 20-85 [iG and highly turbulent both in upstream and downstream, which means 
that the mean free path of accelerated electrons is on the order of their gyro-radius (Bohm limit). This situation 
can be realized both in the age-limited and energy loss-limited model. On the other hand, when the magnetic 
field is nearly perpendicular to the shock normal, which can exist only in the age-limited case, the magnetic field 
is several fiG in the upstream and 14-20 fj,G in the downstream, and the upstream magnetic field is less turbulent 
than the downstream. 
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1. Introduction 

Galactic cosmic rays with energies of less than 10 15 5 eV 
(the "knee" energy) are commonly believed to be gener- 
ated by supernova remnants (SNRs). SN 1006 is one of 
the SNRs which is thought to be an accelerator of such 
high energy particles. Koyama et al. (1995) discovered 
synchrotron X-rays from the rims of this SNR, indicat- 
ing the existence of accelerated electrons with an energy 
more than a few tens of TeV. The detection of TeV 7-rays 
from the northeastern (NE) rim of SN 1006 (Tanimori 
et al. 1998) implies further evidence for the presence of 
high energy particles, since TeV 7-rays arise from the 
Inverse Compton (IC) process, in which cosmic microwave 
background (CMB) photons are up-scattered by high en- 
ergy electrons (Tanimori et al. 2001), or the hadronic pro- 
cess, in which 7r° particles made by collisions between ac- 
celerated and interstellar protons decay into 7-ray pho- 

Send offprint requests to: R. Yamazaki 



tons (Berezhko, Ksenofontov, & Volk 2002; Aharonian & 
Atoyan 1999). 

The mechanism for cosmic ray acceleration has also 
been studied for a long time and the most plausible pro- 
cess is a diffusive shock acceleration (DSA) (Bell 1978; 
Blandford & Ostriker 1978; Drury 1983; Blandford & 
Eicher 1987; Jones & Ellison 1991; Malkov & Drury 2001). 
Many authors have explained the observed properties of 
SN 1006 in the context of the DSA but the conclusions 
are different due to the arbitrary assumptions on unknown 
physical parameters, such as the magnetic field configura- 
tion, the diffusion coefficient, the injection rate, and the 
electron to proton ratio (Ellison, Berezhko, & Baring 2000; 
Berezhko, Ksenofontov, & Volk 2002; Volk, Berezhko, & 
Ksenofontov 2003; Aharonian & Atoyan 1999; Reynolds 
1998; Dyer et al. 2001; Allen, Petre, & Gotthelf 2001; 
Achterberg, Blandford, & Reynold 1998). For example, 
at present the origin of TeV 7-rays can be explained by 
both leptonic and hadronic models. This comes from an in- 
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sufficient theoretical understanding; apart from a globally 
successful picture of the DSA, detailed but important pro- 
cesses, such as the injection or the reflection of accelerated 
particles that determines the above unknown quantities, 
are not well understood. Worse yet, previous observations 
in the hard X-ray band had insufficient spatial resolutions 
to resolve small-scale structures near the shock front, and 
could not strongly constrain the theoretical parameters. 

Recently, Bamba et al. (2003a, b) reported on the re- 
sults for not only spectral but spatial studies of ther- 
mal and non-thermal shock structure in the NE rim of 
SN 1006 with the excellent spatial resolution of Chandra. 
Similar results are also reported by Long et al. (2003) with 
Chandra data. Bamba et al. (2003a) estimated the scale 
length of thermal and non-thermal X-rays in both up- 
stream and downstream of the shock front, which means 
that the direct measurement of the diffusion coefficients 
has become possible. 

In this paper, we show that the important physical 
parameters for the magnetic field in the acceleration site 
can be constrained by the spatial distribution of observed 
nonthermal X-rays. § 2 summarized results of data anal- 
yses by Bamba et al. (2003a). We construct two models 
in § 3; one is assumed that the maximum energy of accel- 
erated electrons are determined by the age of SNR, while 
the other is by the energy loss process such as synchrotron 
or IC cooling. Finally, § 4 is devoted to the discussion on 
validity of our estimation. Throughout the paper, indices 
"u" and "d" represent upstream and downstream, respec- 
tively. 

2. Observed Properties of Nonthermal NE shell 

In this section, we summarize the observed consequence 
of SN 1006, which will be used in the following sections. 

We used the Chandra archival data of the ACIS, which 
has the spatial resolution of 0.5 arcsec on the aim-point, 
on the NE shell of SN 1006 (Observation ID = 00732) 
observed on July 10-11, 2000 with the targeted position 
at (RA, DEC) = (15 h 03 m 51!6, -41 d 51 m 18?8) in J2000 
coordinates. Bamba et al. (2003) and Long et al. (2003) 
report that there are very thin filaments in the hard X- 
ray band image, which must be the acceleration sites of 
high energy electrons. Details of the observation and the 
analysis are found in these papers. 

As shown in Fig. 2 of Bamba et al. (2003), they made 
profiles of six filaments and found that the upstream scale 
length w u ranges between 0.01 and 0.1 pc, while the down- 
stream scale length Wd varies from 0.06 to 0.4 pc using the 
exponential function with adopted distance of 2.18 kpc 
(Winkler et al. 2003). The mean values of w u and are 
0.05 pc and 0.2 pc, respectively. 

Bamba et al. (2003) also fitted the X-ray spectra for 
the six filaments with an sreut model in the XSPEC pack- 
age (Reynolds 1998; Reynolds & Keohanc 1999). The ra- 
dio spectral index of 0.57 was adopted from the result of 
Allen et al. (2001). As a result, the best-fit roll-off fre- 
quency froiioff = 2.6 (1.9-3.3) xlO 17 Hz was derived. The 



quantity z^oiioff is written in terms of magnetic field (B) 
and the maximum energy of accelerated electrons (i? max ) 
as (Reynolds & Keohane 1999) 

^roiioff = 5 x 10 17 Hz ( JL-) ( Emax V . (1) 

Since most of the nonthermal X-ray photons are observed 
in downstream, the synchrotron radiation is mainly due to 
the downstream region. Therefore, it is possible to adopt 
B in Eq. (1) with the downstream magnetic field B d . 

3. Interpretation of the Observed Width of 
Nonthermal X-ray Filaments 

In this section, we explain the observed scale length of 
nonthermal X-ray filaments. Two simple models are con- 
sidered in a context of DSA with a test-particle approx- 
imation. We assume that the spatial distribution of non- 
thermal X-rays coincides with that of the accelerated elec- 
trons to the maximum energy, while thermal X-rays trace 
the spatial profile of a background plasma and hence a 
magnetic field. For a steady state, there is no spatial struc- 
ture of accelerated particles in downstream (Blandford & 
Ostriker 1978). However, one should consider the finite- 
time or energy- loss effect, which makes the spatial profile 
in downstream. 

For simplicity, we assume magnetic fields are spatially 
uniform both in the upstream and downstream at least in 
the nonthermal X-ray emitting region. Since the fraction 
of magnetic pressure to the ram pressure is estimated as 
(B 2 /8ir)/(m H nu^) - 2 x lO^B/lO^G) 2 , where we as- 
sume the number density of thermal plasma n ~1 cm~ 3 
and the shock velocity u s ~ 3 x 10 3 km s _1 , the mag- 
netic pressure does not affect the dynamics of SNR. We 
therefore can adopt the self-similar solution derived by 
Ratkiewicz, Axford, & Mckenzie (1994). Our assumption 
of spatially uniform magnetic field is a good approxima- 
tion in the narrow range around the shock front. 

Since the wide-band spectrum shows a break on the 
X-ray band, electrons accelerated near the maximum en- 
ergy -E max contribute the nonthermal X-ray emission. The 
quantity E max is determined by the age of SNR or by the 
balance of the acceleration and the energy-loss efficiencies. 
To discuss this, we consider three time scales; the acceler- 
ation time scale t acc , the energy loss time scale t\ oss , and 
the age of SN 1006 t age . The energy loss of high energy 
electrons can be neglected when t acc < ti OS s- Then, E max is 
determined by t acc — t agc . On the other hand, E max is lim- 
ited by the energy loss if t acc = t\ oss < t agc . We investigate 
these two cases separately in the following subsections. 

The acceleration time is given by Drury (1983) as 



where u and K are the velocity of bulk flow in the shock 
frame and the diffusion coefficients for accelerated elec- 
trons at the maximum energy, respectively. The diffusion 
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Fig. 1. Contour maps of the functions f(r] u ,8) (left panel) and g(r)d,0) (right). 
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coefficients are given by quasi-linear theory (Skilling 1975; 
Blandford & Eichlcr 1987). Let the mean free path of ac- 
celerated electrons parallel to the magnetic field be a con- 
stant factor J] times the gyro radius r g = E/(cB) in both 
upstream and downstream. Then the diffusion coefficients 
in upstream and downstream can be written in terms of 
?7 U , ?7d, and the angle between upstream magnetic field and 
the shock normal 9 (Jokipii 1987), 



K„ = 



cE„ 



3eB u 



r) u cos 2 9 + 



sin 



1 + 7/5 



(3) 



C.E max . 2 a i 2 • 2 a\— 1 

= - rj d (cos 9 + r sin 9) 



3cB d 



x cos 2 6 + 



r 2 sin 2 i 



(4) 



where r is the compression ratio. In this paper, since we 
assume that the shock is sufficiently strong and that shock 
structure is not affected by the cosmic-ray pressure, r 
should be 4 and the upstream and downstream magnetic 
field are related as 



|l = R(0) : = (cos 2 9 + r 2 sin 2 0)5 . 



(5) 



Upstream and downstream velocities are related as u u = 
ru d = u s , and throughout the paper, we adopt u s = 2.89 x 
10 3 km s _1 (Ghavamian ct al. 2002). It is convenient to 
rewrite diffusion coefficients using Eqs. (3), (4), and (5) as 



Ka 



3cB d 



■f(vu,e) 



3cB d 



9{Vd,0) 



where f(r] u , 9) and 5(r/ u , 9) are given by 
/( % , 9) = ^(cos 2 9 + r 2 sin 2 9) 1 ' 2 ( cos 2 9 + 

(r 2 sin 2 9 
cos 2 9 + - 



(6) 
(7) 

,(8) 
.(9) 



Note that the downstream magnetic field B d is used both 
in Eqs. (6) and (7). Figure 1 shows the contour plots of 
f(r] u ,9) and g(i] d ,9). If 9 ^ 83°, / can be smaller than 
unity in the case of r] u > 1 . The quantity g is smaller than 
0.1 for T] d > 1 and 9 ^ 79°. 

Next, we consider the energy loss time scale. Two 
processes may mainly cause the energy loss of electrons; 
synchrotron radiation and IC effect by cosmic microwave 
background photons. In the following calculations, the lat- 
ter process can be neglected. Then, we can simply write 



Qirm 2 c 
a T EB 2 

1.25 x 10 3 yrs 



^lOOTcVy \10hG) 



(10) 



When one compares t\ oss with i acc , the mean value of the 
magnetic field that accelerated electrons suffer should be 
adopted. We can estimate the mean magnetic field as 



(B 2 ) = aB 2 u + (1 - a)B 2 d 



(11) 



where a = At u /(Ai u + At d ) is the time fraction that 
accelerating electrons are in upstream. For the particles 
in the acceleration process, we can estimate At u / At d ~ 
(K u /u u )/(K d /u d ) = f/(rg) from Eqs. (6) and (7). Then, 
we obtain \ = {R" 2 f +rg) / (f +rg), which ranges between 
and 1, where R is defined in Eq. (5). 

3.1. Age- Limited Case 

First, we investigate the case in which the acceleration 
time is nearly equal to the age of SN 1006 t acc <~ i ag0 . 
This condition implies that the observed nonthermal X- 
rays are emitted by electrons that have been accelerated 
up to now. 

The particles in the diffusive shock acceleration pro- 
cess are transported to the upstream by the diffusion, and 
are advected to the downstream. The diffusion and advec- 
tion time scales to move a scale-length w are t adv — w/u 
and tdif = w 2 /K. In the upstream, the accelerated par- 
ticles can reach up to the place at which the advection 
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is balanced with the diffusion, i.e. t adv = t dif . Therefore, 
the observed upstream- width of nonthermal X-rays can be 
written as 



W 



-f(Vu,0) . 



(12) 



3eB d u s 

In the downstream, if the particles are advected too far 
from the shock front, they cannot return back to the up- 
stream to get further energy. Therefore, the particles being 
accelerated should stay in the region where t d n is smaller 
than t a dv Thus, we obtain the same equation as upstream, 



w d 



K d cE n 



(13) 



3.4 x 10 2 yrs 



tt'u 



(14) 



u d 3eB d u 
Then, Eq. (2) becomes 
t 3r w u + w d ^ 

r — 1 u s \ 0.25 pc 

Since t acc is comparable to the age of the SN 1006, our 
assumption t acc <~ i agc is justified. Eliminating -E max with 
Eq. (1), we can rewrite (12) and (13) as 

- - — - 

^rolloff \ 2 / B d 



0.27 pc 



2.6 x 10 17 Hz 



10 fiG 



f(Vu,0) , 
(15) 



Wd = 1.1 pc 



^rolloff 



2.6 x 10 17 Hz 



Ba 



10 mG 



aim, 



(16) 



respectively. 

Let us take w u ~ 0.05 pc, w d ~ 0.2 pc, and t'roiioff ~ 
2.6 x 10 17 Hz as typical observed quantities (Bamba et al. 
2003a). Then, Eqs. (15) and (16) give B d - 30/ 2 / 3 ~ 
30 ff 2 / 3 fiG. In the case of 6 ~ 0°, the value of B d can be 
determined as follows. We see that since r/ u > 1, / should 
be greater than <~ 1, i.e. B d ^ 30 /xG. Independently of 
this, the condition i acc <, t\ oss gives the upper limit of 
Bd as B d ^ 30 ^G (see Eq. [A.5]). Therefore, we obtain 
B d ~ 30 fiG. This argument has been done in Bamba et al. 
(2003a). On the other hand, as we have mentioned before, 
cases of / ^ 1 can be realized if 9 ^ 80°, i.e. the magnetic 
field is nearly perpendicular to the shock normal. Then 
B d may be smaller than ~ 30 fiG. 

Let us vary observed quantities w u , w d , and ^ ro iioff in 
the range of the observed errors (at 90% confidence level). 
We have six unknown parameters E max , B u , B d , r/ u , r] d , 
and 9. Conditions (1), (5), (12), (13), and t acc ^ ti oss are 
used to relate these quantities. Details are summarized in 
§ A. Figure 2 shows allowed regions of r] u or 7/ d for fixed 
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Fig. 3. Points in the shaded regions surrounded by dotted and dot-dashed lines satisfy Eqs. (1), (5), (17), (18), (19), and 
£ioss < iagc for fixed 8 in the energy loss-limited case. Dotted lines represent Eq. (A. 14) with Wd/w u = 0.6 or 40, while 
dot-dashed lines Eq. (A. 13) with u Iollott = 1.9 x 10 17 or 3.3 x 10 17 Hz. 



9. When we take = 0°, the case of r] u = r/d = 1 (Bohm 
limit in both upstream and downstream) is marginally 
acceptable, since rj should satisfy 1 < rj ^ c/u s ~ 10 2 
(Jokipii 1987). Then the magnetic field has values of 
B u = B d = 20-78 fj,G. The equation (A. 5) shows that 
the maximum value of magnetic fields is achieved when 
observed quantities f ro iioff, w u , and w d have the minimum 
values. On the other hand, when 8 <; 85°, small magnetic 
fields are possible. If we choose w u — 0.1 pc, w d = 0.3 pc, 
and fc'roiioff = 2x 10 17 Hz, then r/ u ~ 10 and jjd ~ 1, the 
magnetic fields are 7?d ~ 4B U ~ 14-20 /zG. This result 
has been suggested by Bamba et al. (2003a). 



The left panel of Fig. 4 represents the allowed region 
of 7? max and B d (see Eqs. [A. 6] and [A.7]). The solid lines 
describe Eq. (1) with i/ roUoff = (1.9-3.3) x 10 17 Hz, while 
the dashed lines a boundary of the region in which a gyro- 
radius of accelerated electrons in the downstream r g = 
E ma ^/(eB d ) oc g~ 1 w d . In order to satisfy 1 < r\ < 10 2 
and all of the other conditions, the quantity g _1 (wd/pc) 
should range between 0.09 and 27. 



3.2. Energy Loss-Limited Case 

If the maximum energy of accelerated electrons E n 
determined by 

tftr.r. ^loss i 



is 



^loss , (17) 

the motion of accelerated particles toward the upstream 
might be obstructed by the energy loss effect as well as the 
advection. Let us consider in the upstream, the energy loss 
time scale icooi 1 in addition to the advection time scale 
tudv — w u /u u and the diffusion time scale idif = w^/K u . 
If idif = ^adv < £ C ooi , the observed width of nonthermal X- 
rays in the upstream is given by w u = K u /u u as well as the 
age- limited case, while in the case of t d n — t coo \ < t &dv , 
w u is given by w u = (K u t coo \)^ . Therefore, we can write 

w u = min{7^ u /u u , (K u t coo i)?} . (18) 

On the other hand, the observed scale length of non- 
thermal X-ray filaments in the downstream w d is deter- 
mined by the cooling time scale as 

w d = max{u d t cooh (K d t coo i)?} . (19) 

1 To avoid the confusion, here we use the notation "t CO oi" 
which is compared with £ adv or t di {, in order to distinguish 
"tioss" used when one compares with i acc or t agc . 
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Fig. 4. The shaded areas indicate the most likely values of magnetic fields just behind the shock front Bd and the maximum 
energy of shock-accelerated electrons E max . The left panel is for the age-limited case, while the right the energy loss-limited 
model. These areas are formed by the two solid lines for a roll-off energy of froiioff = (1.9-3.3) x 10 17 Hz, and the two dashed lines 
from the observed width of nonthermal X-rays; in the age-limited case, a gyro-radius of accelerated electrons in the downstream 
should be r a = £ max /(eB d ) = (0.065-20) x 10~ 2 pc, while in the energy loss-limited case, cooling time of accelerated electrons 
should be t coo i = Wd/wd = (0.81-5.4) x 10 2 yrs. 



We can now use five equations (1), (5), (17), (18), and 
(19) for six unknown parameters E max , B u , Bd, r] u , rjd, 
and 9, and solve these equations with fixed 9 under the 
condition of t\ oss < i ag0 . We summarize a detailed calcu- 
lation procedure in § A. 

Figure 3 shows the results for individual 9. We con- 
sider, as well as the age-limited case, the errors associated 
with the analysis of Bamba ct al. (2003a). When (9 = 0°, 
the case of r? u = rjd = 1 is again acceptable. Then, the 
magnetic field is in the range of B u = Bd = 23-85 /zG. 
If 9 ^ 30°, the downstream magnetic field can be in the 
Bohm limit rjd = 1, then rj u ~ 1-8 and Bd ~ 23-85 fiG. 
However, if 9 is larger than ~ 35°, then rjd <; 10 and 
f] u ^ 10. This implies that the upstream magnetic field is 
more turbulent than the downstream, which seems to be 
unrealistic. 

Indeed, as shown in § A. 2, scale lengths are given by 
w u = K u /u u and Wd — itd*cooi- Then, E max and Bd are 
given by Eqs. (A. 10) and (A. 11), respectively The right 
panel of Fig. 4 represents the allowed region of E max 
and Bd when Wd — 0.06-0.4 pc and v TO noS = (1.9- 
3.3) xlO 17 Hz. All the points in the region satisfy other 
conditions. Vink & Laming (2003) discussed similar argu- 
ments for Cas A. 

4. Discussion 

We have argued the limitations on the model of DSA from 
the recently observed scale length of nonthermal X-ray 
emitting region of NE shell of SN 1006 based on the as- 
sumption that the observed spatial structure of nonther- 
mal X-ray filaments reflects that of accelerated electrons 
to the maximum energy. A test particle approximation 
has been adopted, where the back reactions of acceler- 
ated particles are neglected, although one might have to 
consider corrections due to nonlinear effects in order to ex- 



plain the wide-band spectrum from radio to TeV 7-rays. 
Two models have been discussed; age-limited and energy 
loss-limited models. Note that in each model, the value 
and configuration of a magnetic field in the acceleration 
site can be discussed using only the spatial distribution of 
the synchrotron X-rays, which are emitted by accelerated 
electrons with energy several tens of TeV, and the roll-off 
frequency, which is derived by the fitting of the wide band 
spectrum from the radio to X-rays. 

When a magnetic field is nearly parallel to the shock 
normal, i] u and rjd should be nearly unity, which means 
that the magnetic field is highly turbulent (near the Bohm 
limit), and the magnetic field is in the range of 20-85 jiG. 
Relatively strong parallel magnetic fields are considered by 
several authors (Berezhko et al. 2002; Ellison et al. 2000). 
Derived value of the magnetic field seems to be higher 
than the usual interstellar value of a few micro Gauss. 
However, the mechanism pointed by Lucek & Bell (2000) 
may be able to amplify the magnetic field. These situation 
can be realized both in the age-limited and energy loss- 
limited model. 

When the magnetic field is nearly perpendicular to the 
shock normal, relatively small values of Bd = 4£> u = 14- 
20 /uG are allowed, which can exist only in the age-limited 
model. This value is consistent with the interstellar mag- 
netic field and assumed value in Bamba ct al. (2003a) 
or other papers (Allen ct al. 2001; Reynolds 1998) but 
slightly larger than that derived by CANGAROO obser- 
vation (Tanimori et al. 1998, 2001), which assumes TeV 
7-rays are emitted by the IC process. This discrepancy 
may be solved if the filaments 1-6 are not the main sites 
of TeV 7-rays generated by the leptonic process. The ob- 
served flux of TeV 7-rays which are up-scattered CMB 
photons by synchrotron emitting electrons is written as 

Tic = 2-5 x 10- 2 ( -^-) T synch , (20) 
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where J" S ynch is the observed flux of synchrotron X-rays 
and estimated as 1.8 x 1CT 12 ergs s _1 cm~ 2 in the 0.5- 
10.0 keV band by Bamba et al.(2003a). Then, the contri- 
bution of the IC 7-rays generated in the filaments 1-6 is 
only ~ 0.5% of the whole flux coming from SN 1006. This 
shows that the IC 7-ray emitting sites have smaller mag- 
netic field, hence radiate less synchrotron X-rays than in 
the filaments 1-6. Or the TeV 7-rays might be emitted by 
the hadronic process. In fact the most plausible position of 
the 7-ray emission determined by CANGAROO is located 
north (and outside) of the filament 1 (Tanimori et al. 1998, 
2001). However, the region where the significance is higher 
than half the maximum value extends over ~ 0.2°, which 
is almost the same as the standard deviation of the point- 
spread function of the CANGAROO telescope. We need 
stereoscopic observations of TeV gamma-ray emissions by 
imaging atmospheric Cherenkov telescopes (for example, 
CANGAROO-III, HEGRA, and H.E.S.S.) to establish the 
position and extent of the emission region. 

We briefly discuss how our results may change if 
the test particle approximation is dropped. Although we 
should consider the spatial structure modified, for exam- 
ple, in the upstream precursor region, we simply argue the 
cases of r > 4 that are thought to be caused by nonlin- 
ear effects. Indeed, Berezhko, Ksenofontov, & Volk (2002) 
showed that the present value of total compression ratio 
is about 6. In the age-limited case, we observe that the in- 
ferred magnitude of the upstream magnetic field becomes 
smaller than that of the r — 4 case as r is made larger (up 
to 7). However, the overall shapes and total areas of the 
allowed regions are changed only slightly for any case of 9. 
Thus the nonlinear effect would not be essential. On the 
other hand, in the energy loss-limited case, the allowed 
region in the 77 u -?7d plane becomes narrow and only the 
cases of 9 = 0°-10° can be exist for 4 < r ^ 6 while no 
allowed region for r ~ 7 case. This comes from the fact 
that another restriction, Eq. (A. 24), emerges. The larger 
r, the stronger the constraint because of the less efficient 
advection for fixed shock velocity u s . Furthermore, it can 
be shown that there exists no case in which w,i is given by 
(Kdt coo \)2 , since the condition i acc = t coo \ is incompatible 
with the condition Udt coo i < (i^d^cooi) 3 in the parameter 
range of our interest. 

It is important to determine the magnetic field config- 
urations in order to discuss the acceleration and/or injec- 
tion efficiency. In this paper, using the spatial distribution 
of nonthermal X-rays, we have shown that roughly two 
cases can exist; high and parallel, and low and perpendic- 
ular magnetic field. In the latter case, the back reaction 
of accelerated particles is small and thus a test-particle 
treatment is a good approximation. The magnetic field 
amplification process pointed by Lucek & Bell (2000) does 
not work well. While in the former case, nonlinear effects 
are so efficient that the magnetic field can be large. The 
difference between these cases probably comes from the 
fact that the (proton) injection rate depends strongly on 
the shock obliquity and diminishes as 9 increases (Volk, 



Berezhko, & Ksenofontov 2003). In addition to our result, 
radio polarization data with high spatial resolution may 
become further informations about the magnetic field con- 
figuration. 

In this paper, we have adopted the plane shock ap- 
proximation. For further details, it might be important 
to consider the curvature effect as discussed in Berezhko, 
Ksenofontov, & Volk (2003) in order to produce more re- 
alistic scenario. 
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Appendix A: Detailed calculations 

We summarize a procedure of practical calculations of ?7 U and 
r\ d from observed quantities. Once (/, g) is given for fixed 9, 
the allowed region of (n u ,r)d) is derived by Eqs. (8) and (9). 
Relations of / and g in age-limited case and energy loss-limited 
case are different from each other. 



A.l. Age-limited case 

One can derive from Eqs. (12) and (13) 



9 _ r -i w d 

f W u 



(A.l) 



The condition t acc iioss gives the additional inequalities. 
Substituting Eqs. (6) and (7) into Eq. (2), we derive 



iacc = 4<j>(f + rg) 



cE n 



3e_Bdu| 



(A.2) 



where cj> = <f>(r) — (3/4)r/(r — 1). Energy loss timescale that 
should be compared with t acc is given by 



6nm 2 c 3 f + rg 



(A.3) 



o T E^B\ R->f + rg ' 
Using Eq. (1), we can eliminate -E max and B d , and obtain 

( "* Y ( A 4) 

12.89 x 10 s cm s- 1 J v ' ; 

The quantities / and g should satisfy Eqs. (A.l) and (A.4). 

The condition i acc <, t ioBB gives also the upper limit of B d . 
In the age limited case, one can derive \ = (R~ 2 w u +Wd) / (w u + 
Kid). Then using Eqs. (10) and (14), we derive 



B d <, 30 ^G <j>- 2/3 ( 



frolloff 



-1/3 



2.6 x 10 17 Hz/ 



w u R 2 (6)+w d 
0.25 pc 



-2/3 



,2.89 x 10 s cm s- 1 



2/3 



(A.5) 



where we eliminate E max with Eq. (1). 

Additionally, using Eqs. (1) and (13), one can show E, 
and B d can be expressed in terms of g as 
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Then, t\ oss < i ago reduces 



B d 



31 



9 1 w d 
0.2 pc 

2/3 



41 TeV V3 f-i± 

\4J \2.6x 



rolloff 



1/3 



10 17 Hz, 

u s 



2.89 x 10 s cm s" 



1/3 



(A.6) 



0.2 pc 



froll 



2.6 x 10 17 



off \ 

) 17 Hz J 



1/3 



-2/3 



,2.89 x 10 8 cm s- 1 



-2/3 



(A.7) 



respectively. 



A. 2. Energy loss-limited case 

For our parameters we consider, w u and uid are given by 
uiu = K u /u s , 

Wd = Mdtcool • (A. 9) 

We can calculate _E max and Bd from Eqs. (1) and (A. 9) as 

Rolloff \ 2/3 

10 17 Hz/ 



V4/ \2.6x 
0.2 pc J \ 



Bd = 35 fiG ( - 



0.2 pc 



-2/3 



-2/3 



2.89 x 10 8 cm s- 1 

1/3 



-1/3 



(A.10) 



2.6 x 



^rolloff \ " 

x 10 17 Hz/ 



2.89 x 10 8 cm s 



2/3 



(A.11) 



Then, from Eqs. (6) and (A. 8), we derive 



/ = 4.8^(1 
w d V4 



^rolloff 



2.6 x 10 17 Hz, 

2 



V2.89 x 10 8 cm s" 
The condition £ acc = t\ OBB gives 

R(0)- 2 f + rg = 4.8 ^ 



(A.12) 



10 17 Hz 



G 



,2.89 x 10 8 cm s 
Eliminating t'roiioff from Eqs (A.12) and (A. 13), we derive 



(A.13) 



l + rg/f 
R- 2 + rg/f 



< 3.7 



■Wd 

0.2 pc 

Ms 



/ iago 



(A.16) 



V2.89 x 10 8 cm s" 1 / ^1000 yrs^ 

Using Eq. (A. 14), we rewrite the left hand side of Eq. (A.16) 

as 

Substituting Eq. (A. 17) into Eq. (A.16), we derive 
w d + (1 - R(6)- 2 ) 4> ^) 1 wu 



< 0.74 pc 



©"'( 



2.89 x 10 8 cm s 



1000 yrs 



(A.18) 



(A. 8) This equation is always satisfied since w u and w d range 0.01- 
0.1 pc and 0.06-0.4 pc, respectively. 

Finally in order to validate Eqs. (A. 8) and (A. 9), we con- 
firm the conditions 



K u /u s < (Jfutcool) 2 

Wdicooi > (K d t 

cool J 



(A.19) 
(A.20) 



are always satisfied for the parameters of our interest. Using 
Eqs. (1), (6) and (10), Eq. (A.19) can be rewritten as 



/ < 19.0 



2.89 x 10 8 



frolloff 



2.6 x 10 17 Hz 



This equation, together with Eq. (A.12), reduces 

W u 

— < r , 

w d 



(A.21) 



(A.22) 



which is always satisfied since w u = 0.01-0.1 pc and Wd = 0.06- 
0.4 pc. On the other hand, from Eqs. (1), (7), (10), and (A.20), 
we derive 



g < 1.2 



(j) 



2.89 x 10 s cm s" 1 



^rolloff 



2.6 x 10 17 Hz 



.(A.23) 



Combining this equation and Eq. (A.13), we eliminate f ro iioff 
and u s as 



(A.24) 



3 



(A. 14) 



As long as r < 4, this is always satisfied since 1 — (r/4) 1 (f> = 
(r-4)/(r- 1) < 0. 
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